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Abstract 

We propose a method to estimate the probabihty of new physics 
discovery in future high energy physics experiments. Physics simu- 
lation gives both the average numbers < Ni, > of background and < 
Ng > of signal events. We find that the proper definition of the signifi- 
cance is 5i2 = \/<A^s> + <A''{,> — ■v/<A'';,> in comparison with of- 

ten used sigmficances = 4^ and = , < ^ . 

We propose a method for taking into account systematic uncertainties 
related to nonexact knowledge of background and signal cross sec- 
tions. An account of such systematics is very essential in the search 
for supersymmetry at LHC. We propose a method for estimation of 
exclusion limits on new physics in future experiments. We also esti- 
mate the probability of new physics discovery in future experiments 
taking into account systematical errors. 
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1 Introduction 



One of the common goals in the forthcoming experiments is the search for new 
phenomena. In the forthcoming high energy physics experiments (LHC, TEV22, 
NLC, ...) the main goal is the search for physics beyond the Standard Model 
(supersymmetry, Z'-, VF'-bosons, ...) and the Higgs boson discovery as a final 
confirmation of the Standard Model. In estimation of the discovery potential of 
the future experiments (to be specific in this paper we shall use as an example 
CMS experiment at LHC |j|]) the background cross section is calculated and for 
the given integrated luminosity L the average number of background events is 

< N}j >= (T5 • L. Suppose the existence of a new physics leads to the nonzero 
signal cross section as with the same signature as for the background cross section 
that results in the prediction of the additional average number of signal events 

< Ng >= as ■ L for the integrated luminosity L. 

The total average number of the events is < A'^e^ >=< Ng > + < Niy >= 
{as + ab) ■ L. So, as a result of new physics existence, we expect an excess of the 
average number of events. In real experiments the probability of the realization of 
n events is described by Poisson distribution [§] 

f{n,<n>) = ^^e-<^>. (1) 

Here < n > is the average number of events. 

Remember that the Poisson distribution f{n,< n >) gives Q the probability 
of finding exactly n events in the given interval of (e.g. space and time) when the 
events occur independently of one another at an average rate of < n > per the 
given interval. For the Poisson distribution the variance a"^ equals to < n >. So, 
to estimate the probability of the new physics discovery we have to compare the 
Poisson statistics with < n >=< A';, > and < n >=< > + < Ns >■ Usually, 
high energy physicists use the following "significances" for testing the possibility 
to discover new physics in an experiment: 

(a) "significance" Si < > ■rn, 



(b) "significance" 5*2 



< A^. > 

y/< Ns> + < Nb> 



A conventional claim is that for Si (52) > 5 we shall discover new physics (here, 
of course, the systematical errors are ignored). For ^ A'^s the significances Si 
and 5*2 coincide (the search for Higgs boson through the ^ 77 signature). For 
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the case when Ng ~ N^, Si and S2 differ. Therefore, a natural question arises: 
what is the correct definition for the significance Si, S2 or anything else ? 

It should be noted that there is a crucial difference between "future" exper- 
iment and the "real" experiment. In the "real" experiment the total number of 
events N^y is a given number (already has been measured) and we compare it with 
< Nh > when we test the validity of the standard physics. So, the number of possi- 
ble signal events is determined as Ns = N^v— < > and it is compared with the 
average number of background events < >. The fluctuation of the background 
is a fb = y/Nh, therefore, we come to the Si significance as the measure of the dis- 
tinction from the standard physics. In the conditions of the "future" experiment 
when we want to search for new physics, we know only the average number of 
the background events and the average number of the signal events, so we have to 
compare the Poisson distributions P{n, < Nb >) and P{n, < Nb > + < Ns >) to 
determine the probability to find new physics in the future experiment. 

In this paper we estimate the probability to discover new physics in future 
experiments. We show that the proper determination of the significance is S" = 
\/ < Ns > + < Nb > — y/< Nb >. We suggest a method which takes into account 
systematic errors related to nonexact knowledge of the signal and background cross 
sections. We also propose a method for the estimation of exclusion limits on new 
physics in future experiments. Some of presented results has been published in 
our early paper Q. 

The organization of the paper is the following. In the next Section we give a 
method for the determination of the probability to find new physics in the future 
experiment and calculate the probability to discover new physics for the given 
{< Nb >, < Ns >) numbers of background and signal events under the assumption 
that there are no systematic errors. In Section 3 we estimate the influence of the 
systematics related to nonexact knowledge of the signal and background cross 
sections on the probability to discover new physics in future experiments. In 
Section 4 we describe a method for the estimation of exclusion limits on new 
physics in future experiments. In Section 5 we estimate the probability of new 
physics discovery in future experiments. Section 6 contains concluding remarks. 

2 An analysis of statistical fluctuations 

Suppose that for some future experiment we know the average number of the 
background and signal events < Nb >, < A^s >. As it has been mentioned in the 
Introduction, the probability of realization of n events in an experiment is given 
by the Poisson distribution 
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<n> 



n 



P{n, <n>) 



nl 



(2) 



where < n >=< > for the case of the absence of new physics and < n >=< 
Nh > + < Ng > for the case when new physics exists. So, to determine the 
probabihty to discover new physics in future experiment, we have to compare the 
Poisson distributions with < n > = < Nf^ > (standard physics) and < n >=< Ni^ > 
+ < A^<j > (new physics). 

Consider, at first, the case when < Ni, > ::$> 1, < Ng > » 1. In this case 
the Poisson distributions approach the Gaussian distributions Q 



with /i = cj^ and fj, =< > or fj, =< Nf, > + < Ng >. Here n is a real number. 
Note that for the Poisson distribution the mean equals to the variance. 

The Gaussian distribution describes the probability density to realize n events 
in the future experiment provided the average number of events < n > is a given 
number. In Fig.l we show two Gaussian distributions Pq with < n >=< Nb >= 53 
and < n >=< Nb > + < Ng > = lOA ([|1, Table.13, cut 6). As is clear from Fig.l 
the common area for these two curves (the first curve shows the "standard physics" 
events distribution and the second one gives the "new physics" events distribution) 
is the probability that "new physics" can be described by the "standard physics" . 
In other words, suppose we know for sure that new physics takes place and the 
probability density of the events realization is described by curve II {f2{x) = 
Pg{x, < Nb > + < Ng > , < Nb > +Ng >)). The probability k that the "standard 
physics" (curve I (/i(x) = Pg{x, < Nb >,< Nb >))) can imitate new physics (i.e. 
the probability that we measure "new physics" but we think that it is described 
by the "standard physics") is described by common area of curve I and II. 

Numerically, we find that 



PG{n,^i,a'^) 



aV2TT 



1 




(3) 




With a precision defined by the tails (see Section 5). 
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Figure 1: The probability density functions fi,2{x) = Pg{x, fXi^2,cr^) for 
=< Nb >= 53 and //2 =< > + < >= 104. 



Here ai = y/N^ and a2 = y/Nf, + Ng. The transformation of the distributions to 
standard normal distribution and the exploitation of the equality 

Xq - aj _ Xq - (72 

allows one to find the point xq of the intersection of the curves I and II. 

Let us discuss the meaning of our definition (4). For x < xq = (Ti(J2 we have 
fi{x) > f2{x), i-c. the probability density of the standard physics realization is 
higher than the probability density of new physics realization. Therefore for x < xq 
we do not have any indication in favour of new physics. The probability that the 

f2{x)dx. For X > Xq f2{x) > fi{x) 

-oo 

that gives evidence in favour of new physics existence. However the probability 
of the background events with x > xq is different from zero and is equal to /? = 

/■oo 

/ fi{x)dx. So we have two types of the errors. For x < xq we do not have any 

Jxo 

evidence in favour of new physics (even in this case the probabilty of new physics 
realization is different from zero). For x > xq we have evidence in favour of new 
physics. However for x > xo the fiuctuations of the background can imitate new 
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physics. So the probabihty that standard physics can imitate new physics has two 
components a and /3 and it is equal to k = a + /?. If k equals to 1 new physics 
will never be found in the experiment, if k equals to the first measurement 
with probability one has to answer the question about presence or absence of new 
physics (this case is not realized for Poisson distribution). In other words one can 
say that the area of intersection of the probability density functions of the pure 
background and the background plus signal is the measure of the future experiment 
undiscovery potential. 

As follows from formula (4) the role of the significance S plays 



The signal was assumed to be observable if (1 — e) • 100% upper confidence level 
for the background event rate is equal to (1 — e) • 100% lower confidence level for 
background plus signal (e = 0.01 — 0.05). The corresponding significance is similar 
to our significance S\2- The difference is that in our approach the probability k 
that new physics is described by standard physics is equal to 2e. 

It means that for S\2 = 1,2,3,4,5,6 the probability k is correspondingly 
K = 0.31,0.046,0.0027,6.3 • 10 -^ 5.7 • 10-'^,2.0 • 10"^ in accordance with a 
general picture. As it has been mentioned in the Introduction two definitions 

of the significance are mainly used in the literature: Si = ^ ^ Hi and 

V< Nh > ^ 

52 = < s > rpj^g significance 5i2 is expressed in terms of the 

< Ng > + < Nh > 

S1S2 

significances Si and S2 as S12 = 7; tt- 

01 + ^2 

For < Nb >»< Ns > (the search for Higgs boson through if — > 77 decay 
mode) we find that 



It means that for Si = 5 (according to a common convention the 5a confidence level 
means a new physics discovery) the real significance is 512 = 2.5, that corresponds 
to K = 1.24% (Fig.2). 

For the case A''^ = kN^, S12 = ki2S2, where for k = 0.5,1,4,10 the values 
of ki2 are ki2 = 0.55,0.59,0.69,0.77, correspondingly. For not too high values of 
< Afe > and < A";, > + < A's >, we have to compare the Poisson distributions 
directly. Again for the Poisson distribution P(n, < n >) with the area of definition 
for nonnegative integers we can define P{x,< n >) for real x as 




5i2 ~ 0.5 Si w 0.5 52. 



(6) 
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Figure 2: The dependence of k on number of signal events for cases 5*1 = 
5, ^2 = 5 and S12 = 2.5. 



P(x, < n >) = { . ^ ' (7) 

^ ' ^1 P{[x\, <n>), x>0. ^ ^ 



It is evident that 



r 



P{x, < n >)dx = 1. (8) 

I 

So, the generahzation of the previous determination of k in our case is straight- 
forward, namely, k is nothing but the common area of the curves described by 
P{x, < Nb >) (curve I) and P{x, < Ni, > + < Ng >) (curve II) (see, Fig.3). 
One can find that k = a + P, where 



no 



(< JVt > + < iV. >)" .-«a/.>+<K.>i 



n=0 

l-F{2<Nb>+2<Ns> \2no + 2), 



n! 



/?= E >^\ -<N.> = p^2 <N,> |2no + 2), 

n=nci+l 
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Figure 3: The probability density functions fi,2{x) = -P(x,/ii,2) for fii = 
< Nb > = 1 and =< Nb > + < Ns >= 6. 



^(X^l"-) = n / e'h^-'^dt (see, for example, [|]) and no = [— — 

^2_L( — IjQ LTl{ ± 



Numerical results are presented in Tables 1-6. 

As it follows from these Tables for finite values of < A'^j, > and < Nb > the 
deviation from asymptotic formula (4) is essential. For instance, for Ng = 5, 
Nb = l (5i = 5) K = 14.2%. For N, = Nb = 25 (Si = 5) k = 3.8%, whereas 
asymptotically for ^ 1 we find k = 1.24%. Similar situation takes place for 
Ns ~ Nb. 

3 An account of systematic errors related to 
nonexact knowledge of background and sig- 
nal cross sections 

In the previous section we determined the statistical error k (the probability that 
"new physics" is described by "standard physics"). In this section we investigate 
the influence of the systematical errors related to a nonexact knowledge of the 
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background and signal cross sections on the probability k not to confuse a new 
physics with the old one. 

Denote the Born background and signal cross sections as a'^ and a^. An account 
of one loop corrections leads to —>■ cj^(1 + ^i^) and —>■ <Tg{l + 6is), where 
typically ^if, and are 0(0.5). 

Two loop corrections at present are not known. So, we can assume that the 
uncertainty related with nonexact knowledge of cross sections is around 6if, and 
6is correspondingly. In other words, we assume that the exact cross sections lie 
in the intervals (o"^,cr^(l + 26ib)) and (o"^,(T^(1 + 26is)). The average number of 
background and signal events lie in the intervals 

(<iVO>,<iVO>(l + 2,5ib)) (9) 

and 

{<N^>,<N^>{l + 26is)), (10) 

where < Njj >= ■ L, < >= ■ L. 

To determine the probability that the new physics is described by the old one, 
we again have to compare two Poisson distributions with and without new physics 
but in distinction from Section 2 we have to compare the Poisson distributions 
in which the average numbers lie in some intervals. So, a priori the only thing 
we know is that the average numbers of background and signal events lie in the 
intervals (9) and (10), but we do not know the exact values of < A*";, > and 
< Ng >. To determine the probability that the new physics is described by the 
old one, consider the worst case when we think that new physics is described by 
the minimal number of average events 

< >=< Nj^> + < > . (11) 

Due to the fact that we do not know the exact value of the background cross 
section, consider the worst case when the average number of background events 
is equal to < > (1 + 25ih). So, we have to compare the Poisson distributions 
with < n >=< Nj^ > + < iV° >= 

< Njj > (1 + 26u) + (< > -26u < >) and < n >=< Nj^ > (1 + 26u). 
Using the result of the previous Section, we find that for case < > » 1, < 

> ^ 1 the effective significance is 

^There is a problem to determine systematic uncertainty probability distributions for 
theoretical predictions under consideration. 
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Sus = 7< NS> + < NO > - ^< A^o > {l + 26u). (12) 

For the limiting case 6n, 0, we reproduce formula (5). For not too high 
values of < > and < >, we have to use the results of the previous section 
(Tables 1-6). 

As an example consider the case when 5ib = 0.5, < A'^, >= 100, < AT^ >= 50 
(typical situation for sleptons search). In this case we find that 

S,= , =8.2 

V < iV, > + < A^b > 

Si2 = ^/ <Nb> + <N7>- V <Nb> = 5.2, 
Si2s = V<Nb> + <Ns> - V2 <Nb > = 2.25. 



The difference between CMS adopted significance S2 = 8.2 (that corresponds 
to the probability k = 0.24 • 10""*^^) and the significance S'i2s = 2.25 taking into 
account systematics related to nonexact knowledge of background cross section 
is factor 3.6. The direct comparison of the Poisson distributions with < Nj, > 
(1 + 2(5ife) = 100 and < AT^ > (1 + 2Sii,)+ < Ns,eff > ( < Ns,eff >=< Ns > 
-2(5ib < Nb >= 50) gives Ks = 0.0245. 

Another example is with < A^^ >= 28, < Nb >= 8 and 6ib = 0.5. For such 
example we have = 9.9, ^2 = 4.7, = 3.2, 5125 = 2.0, Ks = 0.045. 

So, we see that an account of the systematics related to nonexact knowledge 
of background cross sections is very essential and it decreases the LHC SUSY 
discovery potential. 



4 Estimation of exclusion limits on new physics 

In this section we generalize the results of the previous sections to obtain exclusion 
limits on signal cross section (new physics). 

Suppose we know the background cross section ab and we want to obtain 
bound on signal cross section ag which depends on some parameters (masses of new 
particles, coupling constants, ...) and describes some new physics beyond standard 
model. Again as in Section 2 we have to compare two Poisson distributions with 
and without new physics. The results of Section 2 are trivially generalized for the 
case of the estimation of exclusion limits on signal cross section and, hence, on 
parameters (masses, coupling constants, ...) of new physics. 
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Consider at first the case when < Nb >= at ■ L » 1, < Ng >= ag ■ L ^ 1 
and the Poisson distributions approach the Gaussian distributions. As it has been 
mentioned in Section 2 the common area of the Gaussian curves with background 
events and with background plus signal events is the probability that " new physics" 
can be described by the "standard physics". For instance, when we require the 
probability that "new physics" can be described by the "standard physics" is less 
or equal 10% (S'12 in formula (5) is larger than 1.64) it means that the formula 



^/<N^^TTTN7> - ^/TN^ <1M (13) 

gives us 90% exclusion limit on the average number of signal events < Ns >. In 
general case when we require the probability that "new physics" can be described 
by the "standard physics" is more or less e the formula 

V< Nb> + < Ns> - V< Nb > < 5(e) (14) 

allows us to obtain 1 — e exclusion limit on signal cross section. Here ^(e) is 
determined by the formula (4) ^ i.e. we suppose that e = k. It should be stressed 
that in fact the requirement that " new physics" with the probability more or equal 
to e can be described by the "standard physics" is our definition of the exclusion 
limit as (1 — e) probability for signal cross section. From the formula (14) we find 
that 




For the case of not large values of < A'^;, > and < Ng > we have to compare the 
Poisson distributions directly and the corresponding method has been formulated 
in Section 2. As an example in Table 7 we give 90% exclusion limits on the signal 
cross section for L = 10'^pb~^ and for different values of background cross sections. 

Formulae (14), (15) do not take into account the influence of the systemati- 
cal errors related to nonexact knowledge of the background cross sections on the 
exclusion limits for signal cross section. To take into account such systematics 
we have to use the results of Section 3. The corresponding generalization of the 
formulae (14) and (15) is straightforward, namely: 



\/< iVs > + < iV, > - J< iVs > (1 + 24n) < S(e), (16) 



,,<^ + ,5(W^5(l±^ + 2*,.. (17) 



^Note that 5(1%) = 2.57, S{2%) = 2.33, 5(5%) = 1.96 and 5(10%) = 1.64 
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Remember that 6ib describes theoretical uncertainty in the calculation of the 
background cross section. As an example, in Table 8 we give 90% exclusion limits 
on the signal cross section for L = W'^pb~^, 25ib = 0.25 and for different values 
of background cross sections. 

Note that in refs.|^, |l^] different and strictly speaking "ad hoc" methods to 
derive exclusion limits in future experiments has been suggested. As is seen from 
Fig. 4 the essential differences in values of the exclusion limits take place. Let us 



compare these methods by the use of the equal probability test [11 1 



„.J.Hernonde2 et ol, 19] oppro 
Typictil experir-ient [1 G] opprc 



Expected background evt 



Figure 4: Estimations of the 90% CL upper limit on the signal in a future 
experiment as a function of the expected background. The method proposed 
in ref. [10] gives the values of exclusion limit close to "Typical experiment" 
approach. 



In order to estimate the various approaches of the exclusion limit determina- 
tion we suppose that new physics exists, i.e. the value < A^s > equals to one of 
the exclusion limits from Fig. 4 and the value < A^ > equals to the corresponding 
value of expected background. Then we apply the equal probability test to find 
critical value no for hypotheses testing in future measurements. Here a zero hy- 
pothesis is the statement that new physics exists and an alternative hypothesis is 
the statement that new physics is absent. After calculation of the Type I error a 
(the probability that the number of observed events will be equal to or less than 
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the critical value no) and the Type II error /? (the probability that the number of 
observed events will be greater than the critical value uq in the case of absence of 
new physics) we can compare the methods. In Table 9 the comparison result is 
shown. As is seen from this Table the "Typical experiment" approach [^0| gives 
too small values of exclusion limit. The difference in the 90% CL definition is 
the main reason of the difference between our result and the exclusion limit from 
ref. 0|. We require that e = k. In ref P] the criterion for determination exclusion 

limits: f3 < A and < e is used, i.e. the experiment will observe with prob- 

l- P 

ability at least 1 — A at most a number of events such that the limit obtained at 
the 1 — e confidence level excludes the corresponding signal 0. 



5 The probability of new physics discovery 

In section 2 we determined the probability k that "new physics" can be described 
by the "standard physics". But it is also very important to determine the prob- 
abilty of new physics discovery in future experiment. According to common defi- 
nition 1^] the new physics discovery corresponds to the case when the probability 
that background can imitate signal is less than 5a or in terms of the probabilty 
less than 5.7 • 10^'' (here of course we neglect any possible systematical errors). 

So we require that the probability /3(A) of the background fluctuations for 
n > n(A) is less than A, namely 

oo 

/3(A) = J2 P{<Nb>,n)<A (18) 

n=reo(A)+l 

The probability 1 — a{A) that the number of signal events will be bigger than 
no (A) is equal to 

oo 

l-a(A)= J2 P{< Nb> + < Ns>,n) (19) 

n=no(A)+l 

It should be stressed that A is a given number and a(A) is a function of A. 
Usually physicists claim the discovery of phenomenon Q if the probability of the 
background fluctuation is less than 5a that corresponds to A^j^ = 5.7 • 10^"^ ^. So 

^If we define e as normalized n {e ~ k ~ ) we have the result close to ref. |Q. For 

example, k — 0.17 corresponds to e = 0.0929, i.e. 1 — e w 0.9. 

^The approximation of Poisson distribution by Gaussian for tails with area close to or 
less than A^is is wrong. 
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from the equation (18) we find rao(A) and estimate the probabilty 1 — a(A) that 
an experiment wih satisfy the discovery criterion. 

As an example consider the search for standard Higgs boson with a mass mh = 
110 GeV at the CMS detector. For total luminosity L = 3 • 10'^pb-'^{2 • W'^pb^'^) 

one can find 0] that < Nb >= 2893(1929), < N, >= 357(238), = f ^ = 

6.6(5.4). Using the formulae (18, 19) for A^j,, = 5.7 • 10^^ (5cr discovery criterion) 
we find that 1 — a(Arfjs) = 0.96(0.73). It means that for total luminosity L = 
3 • 10'*p6~^(2 • 10^p6~^) the CMS experiment will discover at > 5a level standard 
Higgs boson with a mass rrih = 110 GeV with a probabilty 96(73) percent. 

An account of uncertainties related to nonexact knowledge of background cross 
section is straightforward and it is based on the results of Section 3. Suppose 
uncertainty in the calculation of exact cross section is determined by parameter 5, 
i.e. the exact cross section lies in the interval {ab,(Tb{l + 5)) and the exact value 
of average number of events lies in the interval (< A*";, >, < A";, > {1 + 6)). Taking 
into account formulae (18) and (19) we have the formulae 

oo 

/3(A) = Pi<Nb> il + 6),n)<A (20) 

n=no(A)+l 
oo 

l-a(A)= Yl P{< Nb> + < Ns>,n) (21) 

n=no(A)+l 

As an application of formulae (20,21) consider the case < A^^ >=< A";, >= 100 
(typical case for the search for supersymmetry at LHC). For such values of < A^g > 
and < Affe > we have Si = 10, ^2 = 7.1, Su = 4.1. For 6= 0., 0.1, 0.25, 
0.5 we find that 1 - a{Adis) = 0.9998, 0.9938, 0.8793, 0.1696, correspondingly. 
So, we see that the uncertainty in the calculations of background cross section is 
extremely essential for the determination of the LHC discovery potential. Some 
other examples are presented in Tables 10-15. 

Let us consider the random variable "luminosity of 5cr discovery claim" for 
predicted phenomenon in future experiment. Fig. 5 illustrates the behaviour of this 
value for above example < Ns >=< Nb >= 100 at integrated luminosity W^pb~^. 
As follow from Fig. 5(b) we can point out average luminosity of 5(T discovery claim 
L = 0.3287 • lO^pb^^ and estimate the accuracy of this prediction. As seems it is 
very important parameter for comparison of proposals of future experiments. 
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Figure 5: The cumulative distribution function (a) and the behaviour of 
the probabihty distribution (b) of the random variable "luminosity of 5a 
discovery claim" (< Ns >=< Nb >= 100 at integrated luminosity 10^pb~^). 



15 



6 Conclusions 



In this paper we determined the probabihty to discover the new physics in the 
future experiments when the average number of background < Nb > and signal 
events < A^^ > is known. We have found that in this case the role of significance 
plays 5i2 = \/ < Nb > + < Ng > — Nf, > in comparison with often used ex- 
pressions for the significances 5i = ^ ^ ^ = and So = <^ s > 

yj< Nb > y/<Ns> + <A^b> 

For < A'^s > < < iVfe > we have found that = 0.55i = 0.5S'2. For not 

too high values of < A'^s > and < Nb >, when the deviations from the Gaussian 
distributions are essential, our results are presented in Tables 1-6. We proposed a 
method for taking into account systematical errors related to the nonexact knowl- 
edge of background and signal events. An account of such kind of systematics is 
very essential in the search for supersymmetry and leads to an essential decrease 
in the probability to discover the new physics in the future experiments. We also 
proposed methods for the estimation of exclusion limits on new physics and the 
probability of the new physics discovery in future experiments. 

We are indebted to M.Dittmar for very hot discussions and useful questions 
which were one of the motivations to perform this study. We are grateful to 
V.A.Matveev for the interest and useful comments. This work has been supported 
by RFFI grant 99-02-16956. 
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Table 1: The dependence oi k, on < Ng > and < Nb > for Si — 5 





< Nh > 




5 


1 


0.1423 


If) 


4 


n 0828 


15 


9 


0.0564 


20 


16 


0.0448 


25 


25 


0.0383 


30 


36 


0.0333 


35 


49 


0.0303 


40 


64 


0.0278 


45 


81 


0.0260 


50 


100 


0.0245 


55 


121 


0.0234 


60 


144 


0.0224 


65 


169 


0.0216 


70 


196 


0.0209 


75 


225 


0.0203 


80 


256 


0.0198 


85 


289 


0.0193 


90 


324 


0.0189 


95 


361 


0.0185 


100 


400 


0.0182 


150 


900 


0.0162 


500 


10^ 


0.0135 


5000 


10'' 


0.0125 
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Table 2: The dependence of k on < A^^ > and < AT^ > for 5*2 ~ 5. 



<iV, > 


<Nb> 


K 


26 


1 


0.15 


10 


29 


4 


0.14 


10 


33 


9 


0.44 


10 


37 


16 


0.99 


10 


41 


25 


0.17 


10 


45 


36 


0.26 


10 


50 


49 


0.31 


10 


55 


64 


0.36 


10 


100 


300 


0.74 


10 


150 


750 


0.89 


10 



Table 3: < Ng >= -• < >. The dependence of k on < AT^ > and < Nb >. 



<Ns> 


<Nb> 


K 


50 


250 


0.131 






100 


500 


0.033 






150 


750 


0.89 • 


10" 


-2 


200 


1000 


0.25 • 


10" 


-2 


250 


1250 


0.74- 


10" 


-3 


300 


1500 


0.22 • 


10" 


-3 


350 


1750 


0.65 • 


10" 


-4 


400 


2000 


0.20 • 


10" 


-4 
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Table 4: < Ng >= — < >. The dependence oi k on < Ng > and 



<Ns> 


<Nb> 




50 


500 


0.275 


100 


1000 


0.123 


150 


1500 


0.059 


200 


2000 


0.029 


250 


2500 


0.015 


300 


3000 


0.75 • 10^2 


350 


3500 


0.38 • 10-2 


400 


4000 


0.20 • 10-2 


450 


4500 


0.11 • 10-2 


500 


5000 


0.56 • 10-=^ 
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Table 5: < Ng >=< Ni^ >. The dependence of k on < Ng > and < Nh >. 



<Ns> 


<Nb> 


K 


2. 


2. 


0.561 


4. 


4. 


0.406 


6. 


6. 


0.308 


8. 


8. 


0.239 


10. 


10. 


0.188 


12. 


12. 


0.150 


14. 


14. 


0.121 


16. 


16. 


0.098 


18. 


18. 


0.079 


20. 


20. 


0.064 


24. 


24. 


0.042 


28. 


28. 


0.028 


32. 


32. 


0.019 


36. 


36. 


0.013 


40. 


40. 


0.87 • 10-2 


50. 


50. 


0.34 • 10-2 


60. 


60. 


0.13 • 10-2 


70. 


70. 


0.52 • 10-3 


80. 


80. 


0.21 • 10-3 


100. 


100. 


0.33 • 10-^ 
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Table 6: < Ng >= 2- < >. The dependence of «; on < iV^ > and < Nb >. 



<Ns> 


<Nb> 


K 


2. 


1. 


0.463 


4. 


2. 


0.294 


6. 


3. 


0.200 


8. 


4. 


0.141 


10. 


5. 


0.102 


12. 


6. 


0.073 


14. 


7. 


0.052 


16. 


8. 


0.037 


18. 


9. 


0.027 


20. 


10. 


0.020 


24. 


12. 


0.011 


28. 


14. 


0.59 • 10-2 


32. 


16. 


0.33 • 10-2 


36. 


18. 


0.18-10-2 


40. 


20. 


0.10 • 10-2 


50. 


25. 


0.23 • lO--'^ 


60. 


30. 


0.56 • 10-^ 



Table 7: 90% exclusion limits on signal cross section for L = 10^pb~^ and 
for different background cross section (everything in pb). The third column 
gives exclusion limit according to formula (15). 







(Ts (continuous limit) 


10^ 


1.041 


1.038 


102 


0.329 


0.328 


10 


0.104 


0.104 


1 


0.033 


0.033 


0.1 


0.011 


0.011 


0.01 


0.0036 


0.0035 


0.001 


0.0013 


0.0013 


0.0001 


0.00060 


0.00060 
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Table 8: 90% exclusion limits on signal cross section for L = lO'^pb^^, 2Sib = 
0.25 and for different background cross section (everything in pb). The third 
column gives exclusion limit according to formula (17). 



0-6 


0-s 


o"s (continuous limit) 


10^ 


251.25 


251.16 


10^ 


25.37 


25.37 


10 


2.62 


2.62 


1 


0.29 


0.29 


0.1 


0.037 


0.037 


0.01 


0.0064 


0.0064 


0.001 


0.0017 


0.0017 


0.0001 


0.00064 


0.00066 



Table 9: The comparison of the different approaches to determination of the 
exclusion limits. The a and the /3 are the Type I and the Type II errors 
under the equal probability test. The k equals to the sum of a and /?. 







this 


paper 








ref. 


[9] 








ref. 


[10] 






Ns 


a 


P 




K 


Ns 


a 


/3 




K 


Ns 


a 


P 


K 


1 


6.02 


0.08 


0.02 





10 


4.45 


0.09 


0.08 





17 


3.30 


0.20 


0.08 


0.28 


2 


7.25 


0.05 


0.05 





10 


5.50 


0.13 


0.05 





18 


3.90 


0.16 


0.14 


0.30 


3 


8.32 


0.07 


0.03 





10 


6.40 


0.09 


0.08 





18 


4.40 


0.14 


0.18 


0.32 


4 


9.20 


0.05 


0.05 





10 


7.25 


0.13 


0.05 





18 


4.80 


0.23 


0.11 


0.34 


5 


10.06 


0.07 


0.03 





10 


7.90 


0.10 


0.07 





17 


5.20 


0.20 


0.13 


0.34 


6 


10.67 


0.06 


0.04 





10 


8.41 


0.09 


0.08 





18 


5.50 


0.19 


0.15 


0.34 


7 


11.37 


0.05 


0.05 





10 


9.00 


0.08 


0.10 





18 


5.90 


0.17 


0.17 


0.34 


8 


12.02 


0.07 


0.03 





10 


9.70 


0.10 


0.06 





17 


6.10 


0.17 


0.18 


0.35 


9 


12.51 


0.06 


0.04 





10 


10.16 


0.09 


0.07 





17 


6.40 


0.16 


0.20 


0.36 


10 


13.04 


0.05 


0.05 





10 


10.50 


0.09 


0.08 





17 


6.70 


0.22 


0.14 


0.36 


11 


13.62 


0.04 


0.06 





10 


10.80 


0.08 


0.09 





18 


6.90 


0.21 


0.15 


0.36 
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Table 10: The dependence of 1 — a{Adis) on < Ng > and < Nj, > for S'l = 5 
and different values of S. 



<Ns> 


<Nh> 


6 = 0.0 


S = 0.1 


6 = 0.25 


6 = 0.5 


5 


1 


0.0839 


0.0839 


0.0426 


0.0426 


10 


4 


0.1728 


0.1174 


0.0765 


0.0288 


15 


9 


0.2323 


0.1321 


0.0678 


0.0132 


20 


16 


0.2737 


0.1783 


0.0609 


0.0071 


25 


25 


0.3041 


0.1779 


0.0424 


0.0020 


30 


36 


0.3273 


0.1480 


0.0315 


0.0007 


35 


49 


0.3456 


0.1502 


0.0192 


0.0001 


40 


64 


0.3973 


0.1305 


0.0125 


0.00003 


45 


81 


0.4064 


0.1157 


0.0068 


0.000004 


50 


100 


0.4140 


0.1042 


0.0040 




55 


121 


0.4205 


0.0950 


0.0019 




60 


144 


0.4261 


0.0876 


0.0010 




65 


169 


0.4309 


0.0723 


0.0004 




70 


196 


0.4352 


0.0606 


0.0002 




75 


225 


0.4389 


0.0516 


0.0001 




80 


256 


0.4638 


0.0444 


0.00003 




85 


289 


0.4657 


0.0387 


0.00001 




90 


324 


0.4674 


0.0306 


0.000005 




95 


361 


0.4689 


0.0245 


0.000002 




100 


400 


0.4703 


0.0199 






150 


900 


0.5041 


0.0015 
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Table 11: The dependence of 1 — ^(A,^^^) on < iV^ > and < Nj, > for 5*2 ~ 
and different values of 6. 



<N,> 


<Nb> 


5 = 0. 


5 = 0.1 


5 = 0.25 


5 = 0.5 


26 


1 


0.9999 


0.9999 


0.9998 


0.9998 


29 


4 


0.9983 


0.9968 


0.9940 


0.9825 


33 


9 


0.9909 


0.9779 


0.9524 


0.8423 


37 


16 


0.9725 


0.9473 


0.8491 


0.5730 


41 


25 


0.9418 


0.8806 


0.6606 


0.2457 


45 


36 


0.9016 


0.7622 


0.4705 


0.0848 


50 


49 


0.8774 


0.7058 


0.3208 


0.0222 


55 


64 


0.8752 


0.6206 


0.2161 


0.0057 


100 


300 


0.7155 


0.1307 


0.0002 




150 


750 


0.6599 


0.0119 







Table 12: < A^^ >= |- < Nb >. The dependence of 1 - a{Adis) on < K 
and < Nf, > for different values of 5. 



<Ns> 


<Nb> 


5 = 0. 


5 = 0.1 


50 


250 


0.0408 


0.0004 


100 


500 


0.3032 


0.0030 


150 


750 


0.6599 


0.0119 


200 


1000 


0.8905 


0.0301 


250 


1250 


0.9735 


0.0629 


300 


1500 


0.9947 


0.1127 


350 


1750 


0.9992 


0.1767 


400 


2000 


0.9999 


0.2595 
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Table 13: < A^^ >= j^- < Nf, >. The dependence of 1 — a{Adis) on < Ng > 
and < Nh >. 



<Ns> 


<Nh> 


5 = 0. 


50 


500 


0.0043 


100 


1000 


0.0424 


150 


1500 


0.1478 


200 


2000 


0.3223 


250 


2500 


0.5177 


300 


3000 


0.6955 


350 


3500 


0.8270 


400 


4000 


0.9093 
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Table 14: < Ng > = < Nb>. The dependence of 1 - a{Adis) on < TV^ > and 
< Nb > for different values of 5. 



< Ns> 


< Nh> 


6 = 0. 


S = 0.1 


5 = 0.25 


S = 0.5 


2 


2 


0.0003 


0.0003 


0.0001 


0.000005 


4 


4 


0.0016 


0.0007 


0.0003 


0.00003 


6 


6 


0.0061 


0.0030 


0.0007 


0.00006 


8 


8 


0.0131 


0.0041 


0.0011 


0.0001 


10 


10 


0.0218 


0.0081 


0.0027 


0.0002 


12 


12 


0.0467 


0.0206 


0.0050 


0.0003 


14 


14 


0.0589 


0.0283 


0.0080 


0.0004 


16 


16 


0.0956 


0.0512 


0.0116 


0.0007 


18 


18 


0.1401 


0.0609 


0.0156 


0.0007 


20 


20 


0.1903 


0.0925 


0.0200 


0.0012 


24 


24 


0.3005 


0.1402 


0.0395 


0.0017 


28 


28 


0.4122 


0.2280 


0.0656 


0.0031 


32 


32 


0.5166 


0.2821 


0.0969 


0.0050 


36 


36 


0.6089 


0.3773 


0.1323 


0.0073 


40 


40 


0.7268 


0.4703 


0.1704 


0.0101 


50 


50 


0.8762 


0.6688 


0.3216 


0.0181 


60 


60 


0.9572 


0.8309 


0.4397 


0.0332 


70 


70 


0.9865 


0.9206 


0.5784 


0.0612 


80 


80 


0.9960 


0.9648 


0.7205 


0.0850 


100 


100 


0.9998 


0.9938 


0.8793 


0.1696 
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Table 15: < > = 0.5- < >. The dependence of 1 - a{Adis) on 
< Ng > and < Ni, > for different values of 5. 



<Ns> 


< Nh> 


6 = 0. 


S = 0.1 


6 = 0.25 


2 


4 


0.0001 


0.00002 


0.000005 


4 


8 


0.0003 


0.0001 


0.000009 


6 


12 


0.0010 


0.0002 


0.00003 


8 


16 


0.0017 


0.0005 


0.00004 


10 


20 


0.0040 


0.0009 


0.00005 


12 


24 


0.0071 


0.0012 


0.0001 


14 


28 


0.0111 


0.0023 


0.0001 


16 


32 


0.0156 


0.0025 


0.0002 


18 


36 


0.0207 


0.0039 


0.0002 


20 


40 


0.0341 


0.0056 


0.0003 


24 


48 


0.0589 


0.0099 


0.0005 


28 


56 


0.0886 


0.0192 


0.0008 


32 


64 


0.1424 


0.0259 


0.0011 


36 


72 


0.1796 


0.0402 


0.0013 


40 


80 


0.2442 


0.0575 


0.0021 


50 


100 


0.4140 


0.1042 


0.0040 


60 


120 


0.5692 


0.1947 


0.0074 


70 


140 


0.7187 


0.2762 


0.0118 


80 


160 


0.8250 


0.3820 


0.0195 


100 


200 


0.9456 


0.5765 


0.0408 



28 



